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Abstract. Recently, it was suggested that the map-making procedure, which is 
applied to the time-ordered CMB data by the WMAP team, might be flawed by hot 
pixels. This could lead to a bias in the pixels having an angular distance of about 
141° from hot pixels due to the differential measuring process of the satellite WMAP. 
Here, the bias is confirmed, and the temperature two-point correlation function C(i?) 
is reevaluated by excluding the affected pixels. It is shown that the most significant 
effect occurs in C{d) at the largest angles near ■& = 180°. Furthermore, the corrected 
correlation function C(-(?) is applied to the cubic topology of the Universe, and it is 
found that such a multi-connected universe matches the temperature correlation better 
than the ACDM concordance model, provided the cubic length scale is close to L = 4 
measured in units of the Hubble length. 



PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 
1. Introduction 

One of the most important observational inputs of modern cosmology is provided by the 
cosmic microwave background (CMB). The sky maps produced by the WMAP team [1] 
have found wide applications in cosmological parameter extraction and provided several 
surprises. Besides the low power at large scales, which had already been discovered 
by COBE |2], some strange properties emerged concerning the statistical isotropy, i.e. 
an anomalous alignment between the quadrupole and the octopole [31 H] as well as 
an asymmetry of the CMB fluctuations within the two hemispheres with respect to 
the ecliptic plane [5]. Both anomalies are difficult to reconcile within the standard 
inflationary framework. A possible explanation could be related to the uncertainties 
which inevitably arise in the foreground removal procedures [6]. 

Recently, Liu and Li [71 E] have discussed a possible systematics which lies at the 
heart of the map-making procedure. The main point arises from the fact that WMAP 
measures temperature differences at two points on the sphere separated by about 141° 
which is due to the construction of the probe. They demonstrate in [3 E| that the 
pixels, which lie on a scan ring having diameter 141° with a hot pixel at its centre, 
systematically obtain a lower temperature. They show that the corresponding cross- 
correlation function displays a negative correlation at the crucial angle of 141° in the 
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channel 


Ql, Q2 


VI, V2 


Wl, W4 


W2, W3 


diameter 


144.3° 


140.5° 


139.8° 


142.0° 



Table 1. The diameters of the scan rings in dependece on the channels [S]. 

three frequency bands Q, V, and W thus pointing to a systematic distortion occurring 
in all three bands. In order to estimate the magnitude of the induced temperature 
distortion, they [TJ [8] select the 2000 hottest pixels for a detailed analysis and find a 
temperature distortion in the range between 11/iK and 14/iK which in turn translates 
into a 10% difference in the quadrupole and the octopole moments. This systematics is 
not remedied by a prolonged observation time. It should be noted, however, that the 
WMAP team [T] computes the average temperature of pixels, which do not belong to 
the set of hot pixels, by omitting the hot pixels using the processing mask, so that it is 
not clear at which point of the map-making pipeline this systematics arises. 

In this paper, we construct a modified mask in order to check whether the suggested 
systematics really exists. Furthermore, the two-point temperature correlation function 
C {d) is investigated with respect to this systematics. With the temperature correlation 
function, which is based only on those pixels that are not strongly affected by this 
systematics, we carry out a test of a cubic topology of our Universe. 

2. A modified Liu-Li mask 

The hot pixel systematics does not affect all pixels on the corresponding scan rings in 
the same way. To each hot pixel belongs a well defined scan ring, of course, but due to 
the large number of hot pixels there is an equally large number of affected scan rings, 
which intersect each other, in general. This is due to the large size of the scan rings 
having diameters between 135° and 144° depending on the channel , see Table [1] for 
the individual values. Since the suggested systematics leads to a negative bias on a given 
scan ring, the bias should be more pronounced at the intersections of several affected 
scan rings. In order to test this effect, a modified Liu-Li mask is constructed by the 
following procedure. 

We start with the five-year Ql, Q2, VI, V2, Wl, W2, W3, W4 channel sky maps 
which are not foreground reduced and thus contain the hot pixels. These maps as well as 
the KQ75 and KQ85 masks are provided by the WMAP team [Ij on the LAMBDA home 
page lambda.gsfc.nasa.gov. Then all pixels above a given temperature threshold Tthres 
are selected as hot pixels and their corresponding scan rings with the diameters given in 
Table [1] are computed. In this way 8 maps are obtained whose pixel values contain the 
number A^inter of corresponding hot pixels. From these 8 maps we compute 3 average 
maps, one for each of the three bands Q, V, and W with a mean A'inter- If a given pixel 
obtains more than A^inter intersections from affected scan rings, it gets masked out. This 
mask depends then on both thresholds Tthres and A^mter- In a last step all pixels within 
the KQ75 or KQ85 mask are rejected. The remaining pixels should be those being less 



Hot pixel contamination in the CMB correlation function? 



3 



T , , I , , , I , , , I , , , I , , , I , , r 




300 &T 



Figure 1. The distribution of tlic CMB temperature values 5T in /iK is shown for the 
ILC (5yr) map by taking into account only the pixels outside the KQ75 mask. The 
temperature distribution has a standard deviation of 




Figure 2. The degree of the disturbance by hot pixels is shown for the threshold 
^thres — 1 000 /^K in Galactic coordinates. The pixels are represented in dependence 
on the number A^inter of scan rings with respect to the W band belonging to hot pixels 
with 5T > Tthres up to iVinter — 300 (colour online). 

affected by foregrounds and therefore containing the safest cosmological information. 

The value of the temperature threshold Tthres is chosen to be many times larger 
than the typical temperature fluctuation of the cosmological signal. In Figure [1] the 
distribution of the temperature values 6T is shown for the ILC map based on the five- 
year data [1] by taking into account only the pixels outside the KQ75 mask. The 
temperature distribution has a standard deviation of 66/iK. We choose in the following 
three thresholds of lOOOyuK, 2000yuK, and 4000/iK, respectively, which are very large 
compared to the standard deviation of 66yuK of the CMB signal. In Figure [2] we 
show the distribution of the affected pixels for the W band for a hot pixel threshold 
of Tthres = 1 000 /iK, where the Mollweide projection is used in Galactic coordinates. 
The intensity of the pixels depends on the number A'^inter of hits by scan rings of hot 
pixels. The scale is truncated at iVinter = 300. One observes that the most critical pixels 
lie below and above the Galactic plane opposite to the Galactic centre which is in the 
centre of Figure [2l This behaviour is caused by the large diameter of ~ 141° of the scan 
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Figure 3. The mean temperature 5Tmean is shown for the W band foreground-reduced 
map (Syr). The mean value is computed from the pixels inside the modified Liu-Li 
mask for the W band with T^hj-es = 1 000 /iK but outside the KQ75 mask in dependence 
on the threshold A^inter- In Panel a) the value of (^Tmean is computed from those 
pixels which are related to at least -/Venter tiot pixels. The grey band represents the 2a 
uncertainty of the detector noise in the corresponding area. The shaded region denotes 
the 2a uncertainty obtained from 100000 simulations of ACDM models subjected to 
the same procedure. The dotted curves show the 3(t uncertainties. In Panel b) the 
mean temperature (JTintor is computed from pixels corresponding to exactly iVjnter hot 
pixels. 

rings and the fact that many hot pixels above Tthres are near the Galactic centre. The 
modulus of Galactic latitude of the domain with the largest values of A^inter is given by 
180° — 141° = 39° revealed by the dark regions in Figure O 

In order to verify the negative-temperature bias of the pixels belonging to hot pixel 
scan rings, the mean temperature 5Tmean is computed inside the modified Liu-Li mask 
with Tthres = 1 000 /iK but outside the KQ75 mask. The mean (5Tmean depends on the 
chosen value of the threshold number N-^-^ten i- e. pixels are selected if they are related 
to at least A^inter hot pixels. This procedure tests whether the negative bias increases 
with the number A^inter of scan rings. Figure [3^1 reveals that 5Tmean computed for the 
W band foreground-reduced map indeed drops towards more negative values until a 
value of iVinter — 140 is reached. The pixels with iVinter ^ 140 possess a temperature 
bias of the order of —25 /iK. This is the expected behaviour since with increasing iVjnter 
only the most severely affected pixels are taken into account, i. e. those with the largest 
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Figure 4. As in Figure [3] the mean temperatures ST^eim and (JTjntgr are shown but 
for the V band foreground-reduced map (5yr). 

systematics. Whereas panel a) displays ^Tmean computed from those pixels having at 
least A^intor hot pixels, panel b) shows the mean value (5Tinter computed only from those 
pixels having exactly A^inter hot pixels. The number of the latter ones is lower and thus 
the curve possesses larger statistical fluctuations. But it nevertheless demonstrates the 
negative bias for pixels with A^inter above 100. 

In order to check the significance of this result, 100 000 simulations of ACDM 
concordance models are generated for which 5T^ea.n is computed using the same 
procedure as outlined above. These values scatter symmetrically around a zero mean. 
The shaded region in Figure [3|i denotes the 2a deviation from the zero mean. One 
observes that the W band distortion is over a wide range of values of A^inter larger than 
the 2(7 deviation for ACDM models. This demonstrates that there is a systematics for 
these pixels. 

Figures H] and [5] present the result of the same analysis as shown in Figure [3], but 
for the V and Q bands, respectively. In both frequency bands the same negative bias 
of at least 2a is observed when A^mter is sufficiently large, i. e. when enough hot pixels 
are connected to the considered one. For very large values of A^inter the deviation is 
smaller than 2a but very few pixels occur in that range. The values of A^inter cannot 
be directly compared between the different channels since they possess different hot 
pixels lying above Tthres because of the different foregrounds. In the case of the W band 
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Figure 5. As in Figure [3] the mean temperatures ST^eim and (JTjntgr are shown but 
for the Q band foreground-reduced map (Syr). 

the dust contribution dominates whereas the synchrotron and free-free emission is more 
important in the case of the V and Q bands. 

None of our 100 000 simulations of ACDM concordance models yields a 5Tmean curve 
which always has values below the WMAP result shown in the panel a) of Figures [3] - 
El The same result is obtained for (5Tinter shown in the corresponding panel b). This 
emphasises the peculiarity of the negative-temperature bias. 

3. The correlation hole at ^ = 180° 

Let us now come to the temperature two-point correlation function C{d), which is 
defined as 

C{d) := {5T{n)5T{n')) with n-h' = cos^ . (1) 

It was computed by [2] from the COBE 4yr data and was found to display unexpectedly 
low power at large angles i) > 60°. At these large angles (^(-i?) is close to zero. Only 
near i) = 180° it shows a slight negative power, which is much more pronounced in 
the WMAP five-year data [10]. This anticorrelation of antipodal points is the so-called 
"correlation hole". Copi et al. [10] show the strong dependence of the correlation hole 
on the chosen mask. The correlation function C{'d) possesses the largest antipodal 
anticorrelation if one uses no mask at all, whereas applying the KQ75 mask reduces 
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Figure 6. The correlation function C^d) is computed from the W band foreground- 
reduced map (5yr) using the pixels outside the KQ75 mask. In addition the pixels 
inside the modified Liu-Li mask are rejected for 4 different values of the threshold 

-Winter- 

the anticorrelation, but it is nevertheless present. They conclude that the full-sky 
results seem inconsistent with the cut-sky results. Furthermore, most of the large-angle 
correlation in reconstructed sky maps lies inside the part of the sky that is contaminated 
by the Galaxy. 

We do not take into account the cosmic variance which arises on the theoretical 
side since the cosmic initial conditions are assumed only on statistical grounds. Here, 
we are interested in the correlation function C{'d) observed from our special observer 
point within our Universe. 

We now show that the remaining negative correlation near d = 180° is caused by the 
pixels lying within the modified Liu-Li mask. To that aim we apply our modified Liu- Li 
mask with Tthres = 1 000 fiK together with the KQ75 mask to the W band foreground- 
reduced map (Syr) for several choices of the intersection threshold A^inter- The result 
shown in Figure [6] reveals a pronounced correlation hole for large values of A'^inter, i- e. 
for those cases where only a small number of further pixels is rejected after applying 
the KQ75 mask. The curve belonging to A^inter = 200 displays almost the usual result 
obtained by using only the KQ75 mask, whereas the one belonging to A^inter = 60 shows 
even positive correlations near i) = 180°. The four cases presented in Figure [6] also 
demonstrate that the correlation is at other angular separations nearly independent 
on the value A'inter except for the interval [30°, 50°]. However, the strongest sensitivity 
occurs at "i? = 180°. 

Let us now address the question how probable such a behaviour is with respect to 
the ACDM concordance model. In order to compute the probability that C{^) remains 
nearly unchanged for 50° < '& < 165°, define the difference 

ACi^) := CkQ75+Ll(^9) - C7kQ75(^9) (2) 

between C('(9) using the KQ75 mask and the modified Liu-Li mask for A^inter = 60, which 
is denoted as Ckq75-i-ll ("i?) , and C{^) using only the KQ75 mask denoted as Ckq75(^?)- 
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Figure 7. The weighted temperature correlation difference I{L) is computed from 
the four correlation functions C{'d) shown in Figure [6l The length L denotes the size 
of the cubic topology measured in units of the Hubble length. 

The integrated quantity 

D := dcos^ lAC'WI (3) 

J-d=50° 

measures how strongly C{i)) is affected by the suspected pixels. The value obtained 
from the W band data is surprisingly low, indeed from our 100 000 ACDM concordance 
models only 495 models possess a smaller value, i.e. -Dacdm < -Dwmap- Thus, 99.5% 
of the ACDM models show a stronger change in C{i)) than the considered WMAP 
data. Concerning the correlation hole, only 3111 models satisfy ACacdm(180°) > 
ACwMAp(180°), i.e. the large increase of C('(9) at = 180° is also unusual at the 3 
percent level. If one requires both conditions simultaneously not a single ACDM model 
remains. This emphasises the unusual behaviour of the suspected pixels. 

4. The cubic topology for the modified correlation function 

The surprisingly low CMB large-angle power of the correlation function C {i}) , which was 
discovered in |2j, is at variance with the ACDM concordance model as has been found 
in yj] and recently emphasised in [12] and [10]. One of the possible explanations for 
this unexpected behaviour is provided by models of the Universe whose spatial section is 
represented by a multi-connected space form. For a general introduction to the topology 
with respect to cosmology, see e.g. [13], HH |15l |16l |171 [18]. In these models exists 
naturally a largest wavelength at which perturbations can occur as seeds for structure 
formation. This wavelength is mainly determined by the size of the topological space 
form. The correlation function C('(?) has been studied for several multi-connected space 
forms, and a low large-angle power has been shown to exist in the case of the hyperbolic 
Picard topology [13 [20], the spherical Poincare dodecahedron [21], and the flat toroidal 
universe [T2] . 

In [12] we have analysed the WMAP three-year data with respect to the question 
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whether a flat toroidal universe, i. e. a cubic topology, is compatible with the CMB 
observations. To that aim the correlation function C™°*^'^' (i?) obtained from the cubic 
topology with side length L is compared to the observed correlation function C°^^{'d) 
by using the integrated weighted temperature correlation difference 



It is found in [12] that the results depend on the chosen mask that is applied to the 
CMB data. Applying the kpO mask to the data shows that cubic topologies around 
L = 3.86 are preferred whereas using no mask at all leads to models having a size 
around L = 4.35. {L is given in units of the Hubble length Lh = c/i^o — 4.26Gpc for 
h = 0.704.) This demonstrates the sensitivity on the pixels that are selected for the 
computation of the correlation function C{i}). 

Instead of using the three-year data as in \J2\, we now compute the measure I{L), 
eq. (jl]), using the updated five-year data, where we use the W band foreground reduced 
map. Furthermore, we apply the KQ75 mask together with the four modified Liu-Li 
masks discussed above. The result shown in Figure [7] is that cubic topologies slightly 
smaller than L = 4 lead to the smallest values of I{L), i. e. lead to the best agreement 
with the data. The corresponding values for the ACDM concordance model having 
infinite volume can also be read off from Figure [71 Since the distance to the surface of 
last scattering (SLS) is Lsls = '^V^u — 14.2Gpc where Arj = r/o ~ ''7sls = 3.329 {rj is 
the conformal time), topologies above L 6.6 are so large that the SLS fits completely 
inside the toroidal universe. For larger models the value of I{L) saturates at the infinite 
limit L oo corresponding to the concordance model. For all four values of A^inter 
the correlation function C{'d) is in much better agreement with a cubic universe with 
L slightly below 4 than with the infinite ACDM concordance model. The application 
of the modified Liu-Li mask has little infiuence on the ratio of /(L)//(ACDM) for the 
minimum near L = 4 which is due to the fact that (^(t?) is practically unchanged over the 
large interval from ^ ~ 50° up to -i? ~ 165° and below ^ 30°. Excluding more pixels, 
i. e. decreasing A^inter form 200 to 60, produces a more pronounced second minimum 
around L ^ 2 which lies for A^intor = 60 below the value of the concordance model but 
above the absolute minimum near L = 4. This new feature is caused by the diminished 
power of C{i)) for angles in the interval [30°, 50°] leading to a C{i)) with a power 
suppression already above {} ~ 30° instead of the usual 60°. 

5. Summary 

The presently best full-sky maps of the cosmic microwave background radiation are 
obtained by WMAP whose receivers are differential radiometers measuring the difference 
between two telescope beams separated by about 141°. Liu and Li [?[ [8] discuss the 
possibility that the inclusion of hot pixels in the map-making process causes a bias in 
those pixels corresponding to a scan ring of about 141°. Due to the large diameter of the 
scan rings, the affected pixels lie in a region far from the galactic plane which is usually 
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considered to be scarcely contaminated. Here we confirm the bias found in [TJ E] and, 
furtliermore, compute tlie temperature two-point correlation function C{'d) by taking 
this bias into account. The most dramatic effect of these biased pixels is the so-called 
correlation hole dX d = 180° in C{'d), which is absent by eliminating sufficiently many 
affected pixels. Thus, a correlation function C{'d) is obtained which displays almost no 
correlations above all angular scales larger than 60°. 

In Figures [iHZl we have presented our results for the modified Liu-Li mask using 
the threshold temperature Tthres = 1 000 /xK and employing the KQ75 mask. We have 
checked that the conclusions of this paper remain the same if the KQ85 mask is used 
instead of the KQ75 mask and, furthermore, if the calculations are performed for the 
thresholds 2000/iK and 4000yuK, respectively. 

Using the correlation function C{d) based on the most reliable data, a test for a 
universe with a cubic topology is carried out. It is found that such a model fits C{'d) 
better than the ACDM concordance model if the side length L of the cubic fundamental 
cell is slightly smaller than L = 4 measured in units of the Hubble length. 
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